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Q ' Abstract 

, We examine the calculation of density perturbation at large scales in the inflationary 

' scenario. The formula for its magnitude is reviewed from first principles and applied to 
the original model of extended inflation. Our estimate is an order of magnitude bigger 

. than the earlier ones due to the difference in the time at which the primordial fluctuation 

I is evaluated and to the inclusion of terms neglected earlier in the standard formula for 

'sj- ■ density perturbation. 
(N 

^ 1 Introduction 

Qh, One of the remarkable features of inflationary models is that, besides solving the classical 
Q^l problems of large scale homogeneity and flatness, it also addresses the problem of creating 
the seeds of density inhomogeneity in the early universe. These are the quantum fluctuations 
in energy density during the inflationary epoch, which might have evolved to produce the 
observed density inhomogeneity in the universe. Much work has naturally been devoted to a 
^ i study of the evolution of the initial density perturbation |0, ||, 

I The difficulty here is that of constructing a realistic model ||^. In the original, the 'old' 

inflationary model 0, the bubbles of true phase, nucleated during the strongly flrst order phase 
transition, cannot percolate. The slow-rollover models require too small a coupling constant 
for reproducing the observed density fluctuation to cause problem of reheating. Finally the 
so-called 'chaotic' models require initial conditions having no resemblance to the thermal 
equilibrium condition (or deflnite departures therefrom), usually believed to prevail in the 
early universe. 

The extended inflationary models 0] are similar to the 'old' model with the Einstein 



gravity replaced by the theory of Jordan and of Brans and Dicke . The idea is to take 
advantage of the time dependence of gravitational 'constant' to solve the bubble nucleation 
problem. Though a realistic model still eludes construction [|12| , further work in this direction 



is expected to lead to such a model without violating astrophysical observations p!3| . 

In this work we examine the evaluation of density perturbation in the inflationary epoch 
||14|| . We flrst relate the density perturbation to the quantum fluctuations in the energy density 
operator from flrst principles O, |16| . We then consider the original model of extended 
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inflation in the so-called Einstein frame [|T7[. By restricting the energy density operator to 



terms linear in the scalar field, the quantum fluctuations in the density reduce essentially to 
the two point function of the scalar field. The density perturbation turns out to be an order 
of magnitude bigger than the earlier evaluations [|1^, 0. 

In sec. 2 we derive the formula for the density perturbation from first principles. In sec. 3 
we consider the original model of extended inflation and derive the two point function for the 
scalar field. This is used to evaluate the density perturbation in sec. 4. In sec. 5 we identify 
the points of difference between the present and the earlier evaluations leading to an order of 
magnitude difference. Finally in sec. 6 we summarise our work emphasising these differences. 



2 Density fluctuation formula 

We begin by reviewing the definitions of density inhomogeneity and quantum fluctuation to 



relate them from first principles |15, 16|. Let p{x, t) be the energy density field of the Universe, 
which we assume for simplicity to be confined within a (large) volume V . The mean square 
fluctuation in p(x, t) is defined as 



'5p\ _l( p{x,t)--p{t) V'^ 



(2.1) 



where (■ ■ denotes average over all space and p(t) = {p{x,t))x, the averaged, homogeneous 
background density. One writes 

p{x,t)=p{t){l + 6{x,t)), (2.2) 
where the density contrast 6{x,t) has the Fourier expansion 

5ix,t) = ^j:Sk{t)e-^''- (2-3) 
Then the mean square fluctuation (2.1) becomes 

The so-called fluctuation power per logarithmic interval in wavenumber is defined as 

fSpV P|4(t)P 



(2.4) 



Ph 27r2 



(2.5) 



For our purpose it is useful to consider a related quantity, namely, the mean square mass 
fluctuation on a given length scale. It measures the mass fluctuation within a certain volume 
V by averaging the squared excess mass in it over all points Xq throughout the volume V of 
the universe. To avoid sharp boundary, one smears v with a Gaussian window function. The 
mass within such a smeared sphere placed at xq is 

mi{xo,t) = [ £ye-y"^^'"pixo + y,t). (2.6) 
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The mean square mass fluctuation on length scale / is then 



'^Y _ I ( mi{xQ,t) - mi{t) \\ 



Here mi{t) =< mi{xo,t) >x'o, the average mass obtained by replacing p by p in (2.6). The 
subscript c stands for classical and indicates the phenomenological nature of the evaluation. 
Inserting (2.3) in (2.6) we get 

Now suppose that this density inhomogeneity has a quantum origin in the inflationary 
epoch. Then it should be calculable by considering the quantum fluctuations in the energy 
density operator, Too. This operator will be written explicitly for the specific Lagrangian of 
extended inflation considered in the next section. Here we write it as a sum of two parts, 

Tooix,t)=p{t) + U{x,t), (2.9) 

where p, the classical part, is to be identified with the background density considered above 
and U, the quantum part, depends on the quantum field operator. 
In analogy with (2.6) let us define the mass operator. 



mi 



{t) = J d'xe-^"/^'Xo{x,t), (2.10) 



where the hat on m emphasises that it is an operator. The mean squared fluctuation in mass 
in the smeared sphere is 

fSmV ^ If miit) - mi{t) V\ ^^^^^ 



^ /l,q \\ "^'W 



where (■ ■ ■) denotes the quantum mechanical expectation value in the state with energy den- 
sity p{t), which corresponds to the vacuum of the conventional quantum field theory. The 
subscript q distinguishes it from the phenomenological evaluation in (2.7-8). Inserting the 
Fourier transform of the two point function of U 



{U{x,t)U{x',t)) = / -0^s^-''-^'-'^'^D,{t), (2.12) 



in eqn (2.11), we get 



(5m\^ r d?k _j^2i2Dk{t) 



mj^^^ J (27r)3 p^{t)- 



(2.13) 



If we now identify the spatial averaging in (2.7) with the quantum mechanical expectation 
value in (2.11), we get the desired relation, 

^ . (2.14) 
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allowing us to calculate the density fluctuation during inflation, when it is still within the 
Hubble length (causal horizon). Its evolution as superhorizon-sized perturbation, from the 
time th when it leaves the Hubble radius until its re-entry within it later, is described by the 
constancy of the quantity ( = Skit) /{I +p/p) Assuming radiation dominance at the time 
of re-entry, we flnally get the density perturbation at that time as 




^ (2.15) 



3 Extended inflation in Einstein frame 

Because of the nonminimal coupling of the scalar fleld in the original Brans-Dicke (BD) Action 
|ll| , it is not possible to do calculations with it in the canonical framework of quantum fleld 



theory. It has, however, been shown [|T^ that an appropriate Weyl rescaling can transform 



'2cj + 3 

-nip, 



this action to a form where both the gravity and the kinetic term in the scaled BD fleld are 
canonical. In this so-called Einstein frame, the rescaled action becomes 

S = J d'x^ + \g^''d,^d,^ - (3.1) 

where \l/(x) is the BD fleld in the Einstein frame and 

V lOTT 

nip = being the present value of the Planck mass. This form for the potential function 

results from assuming the matter Lagrangian in the original (Jordan) frame to be dominated 
by the false vacuum energy density, M^. The dimensionless parameter a; appears in the 
original BD Lagrangian. 

We decompose the BD fleld \1/ into a homogeneous, classical fleld if) and a quantum fleld 

^(x,t) = ^/'(t) + 0(f,t). (3.2) 

-> 2 

In the spatially flat FRW metric, ds^ = dt^ — a{tYdx , the classical equations of motion for 
il){t) and the scale factor a{t) are 

a dV 

^ + 3- + — = (3.3) 

a dip 



d\2 Stt /I 



a J ?)m?p \2 



r + V{iP) (3.4) 



In the Einstein frame, the extended inflation resembles a slow-rollover inflation off an expo- 



nential potential |18]. 



The solution to the coupled equations are 0, |T^ 

V^(t) =V'(0)+V^oln(l + Ct), (3.5) 
a(t) = a(0)(l + Ct)P, p = (2cj + 3)/4. (3.6) 
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where 



C = IMLe-m/^o , _ . / (6.; + 5) (2^ + 3) 



qujuip V 327r 

We note that the Hubble length in this model is, 

We now turn to the quantum theory. Expanding the Lagrangian around the classical field 
one gets an (infinite) series in powers of the quantum field 0, though the coefficient of ex- 
pansion (M/mp) is small. Such a Lagrangian is, of course, not perturbatively renormalisable. 
Here we simply reject all higher order terms retaining, in fact, only the lowest (second) order 
terms to get 

S, = J £xdta\t) (^02 - ^(V0)^ - \^?{t)) , (3.8) 

where 



We calculate the Feynman propagator for in the classical background field ip{t). Its 
Fourier transform in 3-space, 

(T0(x,t)0(x'',t')) = / ^e-^^ (^^"-')G,(t,t'), (3.9) 
satisfies the inhomogeneous equation 

The Green's function can be constructed by the familiar procedure. Let /'•^-'(t), f^~\t) be two 
linearly independent mode functions satisfying the homogeneous part of the above equation 
and incorporating the boundary conditions that Z*-^^ contain respectively positive and negative 
frequencies during the initial period for not too low physical momenta |^ . Their normalisation 



is fixed by the value of the Wronskian of the two solutions derivable from (3.10). Then we 
have 

Gfc(t, t') = f+\t)f~\t')e{t - 1') + f^-\t)f^+\t')e{t' - 1). (3.11) 

The equation for the mode functions can be solved exactly. With a change of variable P^ , 
r = (1 + Ctf-v/i^p - 1), it becomes 

(d" 2p 1 d 2(3p-l) \ (3 12) 

It is now easy to cast it in the standard form of Bessel's equation. We get 

f(+'-)(t) - 1 ^ll^lfM (3 13) 



where Hj)-''^^ (z) are the Hankel functions of the first and second kind respectively p5| , ^ and 



/3(3p-l)(p-3) p k 

V = ; ^ , Z = KT 



2{p-l) ' p-laH' 

Two point functions involving and are immediately obtainable from (3.11) and (3.13). 

We also write here the expressions for the energy density and pressure, which can be 
obtained directly from the energy-momentum tensor, 

T,. = d,^d,^ - g^A^g^^'do^'^dp^ - ^(^)), (3-14) 

for the action (3.1). We have already split the energy density as the sum of a homogeneous, 
classical part and a quantum part in (2.9). Inserting (3.2) in (3.14), they are 

p = 1^2 ^ ^ ^mlH\ (3.15) 



U = - — + ^0 



"^^^ =[(6cj + 5)0+ (2cj- l)r-^0]. (3.16) 



7r(2a; + 3)3 dr 

In U we have only retained terms linear in 0. Higher order terms like 0^ would give rise to 
loops and hence divergence in the expression for the density fluctuation. As already stated, 
such terms are multiplied by small coefficients. Thus in any reasonable renormalisation scheme 
, such terms are expected to produce small contributions. We also note the expression for the 
homogeneous pressure, 

P = - = -^^^^mlH\ (3.17) 

4 Evaluation of 6p/ p 

Since the energy density operator f/(0) is linear in 0, it is simple to express the two point 
function of ?7(0) in terms of the mode functions f^^\t). Noting the symmetry of Hankel 
functions, H^^^z) = Hj)^^*(z) for real z, we can write its Fourier transform (2.12) as 

Mt) = -^S^|(6c. + 5)/W(r) + (2c.-l)r-^/W|2, 
n[2uj + 3)3 dr 

^^W^Vwr. (4.1) 



4(2tu- l)(2w + 3)2 \ a{t) 
where 

F{t) = ^{6u + 5)HW(z) + {2u - l)z^Hi'\z) 



Let Xphys{t) be the physical wavelength characterising the density perturbation at time t, 
belonging to the comoving wavenumber k, Xphysif) = 27ia(t)/k. At time th during inflation, 
when it equals the Hubble length, its magnitude is Xphysith) = H^^{th), so that k/{aH)\t^ = 
27T. Then eq.(2.15) becomes 

\ P J H 3v2u; -1 nip 

At time the argument of the Hankel function becomes z{th) = 27rp/ (p — 1), a value large 
enough to justify the use of asymptotic expansion [|^. Retaining terms in z upto the first 
nonleading one, it is simple to evaluate F(th) giving 



F{th) = 2^J{2uj-l){2uj + 3)K{uj), (4.3) 

where K is of order unity. 



^2. . , 3(20c<;2 + 48cj + 27) 



47r2(2cj + 3)2 

To determine th, we follow earlier authors to make certain simplifying assumptions. 
It is assumed that the extended inflation ends at time and instantly gives rise to the radiation 
dominated era with an initial temprerature T ~ M. Further, since the gravitational 'constant' 
ceases to vary in the radiation era, the (classical) BD field ip is set equal to zero from the time 
te onwards |2^. We then get from (3.5), 

M quj nip 



H{te) 2p M ' 

The corresponding wavelengths at time th and at the present time tp can be related as 

Aphys[th)- ^^^^^ aitp)"''''^ ''- 
The scale factors can now be evaluated to give 

M (M/H{th)Y 



(4.4) 



H{th) \M/H{te)J 



TpXphysitp), (4.5) 



where Tp = 2.7 K, the present temperature of the background radiation. Solving for M{th 
and using (4.4) we get 



2p M 



H{th) = M (^^^j {Tp\phys{tp))^ . (4.6) 
Inserting this value in (4.2) we finally get 

=|-irv^[(2cu + 3)/2]5I^(go;)-— (M/mp)^(TpA(g) — . (4.7) 



This result differs from what one finds in the literature 0, |T^ . In particular it is larger than 
the result by Guth and Jain []T^ by a factor of SnK/S ~ 10. The sources of this difference are 
discussed in the next section. 
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5 Comparison with earlier works 



There are two major differences witli earlier works, wliicfi contribute to tlie enliancement of 
our result. One is the time at which the mode functions are evaluated. As already stated, 
the argument of the Hankel functions around the time of Hubble length crossing is sufficiently 
large for the asymptotic expansion to apply. 




{z large). (5.1) 

Then the fluctuation in these oscillatory Fourier components of the quantum field is given by 

k^\f^\t)\^ _( k hV 



27r2 I aH 27r 



(5.2) 



which reduces to just at the time of horizon exit, independently of the value of v. 

But in the literature, the evaluation at the Hubble length crossing is actually carried out 
when the wavelength of the perturbation is much larger than the Hubble length and the 
corresponding mode ceases to oscillate. In this 'frozen' state, the argument of the Hankel 
functions is small enough for expansion around the origin to apply, 

i/^^'2)('2) ^±^^(-)"', {z small). (5.3) 

TT 2 

Further v is assumed very close to | . Under these conditions the fluctuation in the quantum 
field is given by 

The superscript (o) indicates evaluation following literature. Thus our expression (5.2), though 
not of this form in general, does reduce to this form at the time of Hubble length crossing but 
is larger by a factor of (27rp/ (p — 1))^. 

The other difference is the omission of the term in the energy density operator. Thus in 
the literature the energy density fluctuation is related to the field fluctuation by 

Comparing with (4.1) we see that the neglected term is of the same order of magnitude as the 

one retained for z = 27rp/ (p — 1). 

At this point it is simple to obtain the standard expression for the density fluctuation and 
verify the enhancement factor. For the slow-rollover scenario, (dV/dijj) may be estimated from 
the classical equation of motion (3.3) as (dV/dijj) — —3Hip. Then inserting (5.4-5) into (2.15) 
we recover the standard formula, 

' (5.6) 



P 



H 
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on ignoring a factor of 2(p — l)/7rp. It can be immediately evaluated to give 

{'-^Y = .jM^)^ (5.7) 
\p)h "^p 

This standard result is clearly smaller than our result (4.2-3) by the factor St^K/?,. 

6 Discussion 

We review from first principles the formula for the density perturbation due to quantum 
fluctuations in the energy operator. We then evaluate the resulting formula in the context of 
the original model of inflation in the Einstein frame, using the simplifying assumptions made 
by the earlier authors. 

The present evaluation of the density fluctuation differs from the earlier ones in two re- 
spects. One is the time at which the perturbation in the inflationary epoch is evaluated. We 
evaluate it exactly at the time when the characteristic wavelength starts growing bigger than 
the Hubble length. Then the mode functions are still oscillatory and physical quantities admit 
conventional interpretation. In the literature, however, the perturbation is actually evaluated 
when the wavelength has grown several times bigger than the the Hubble length and the mode 
functions have ceased to oscillate. 

Although evaluated at two different times with qualitatively different mode functions, the 
density fluctuations turn out to have the same, nearly scale independent, spectum. It is 
interesting to note that our evaluation leaves the index v of the Hankcl function free, while 
that in the literature requires v to be equal to | to arrive at this spectrum. But the two 
evaluations do make a difference in the magnitude of the fluctuation, our result being greater 
by a factor of 27r. 

The other point of difference with the earlier evaluation is the inclusion of the term in 
Too. In our way of evaluation, this term is as important as the term. These two sources of 
difference go to make our result larger by about a factor of 10. 

We hope to have made it clear that these sources of enhancement in our calculation over 
the earlier ones, although discussed here in the context of a specific model, is, in fact, quite 
general. But to arrive at the physically correct magnitude of {5p/p)H, we need a better 
knowledge of the evolution of the density perturbation than what the approximate constancy 
of ( suggests. Pending this investigation, a comparison of our result with the earlier ones 
indicates the magnitude of uncertainty in such a calculation. 

One of us (SM) would like to thank Professor Paul J. Stcinhardt for earlier discussions on 
related topics. He also gratefully acknowledges the hospitality at the Institute for theoretical 
Physics, University of Bern, which made these discussions possible. 
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